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CONTTROLLABILITY  AND  STABILIZABILITY  OF 
REGULAR  SINGULAR  LINEAR  SYSTEMS  WITH  CCWSTANT  COEFFICIENTS 

by 

C.  E.  Langenhop 

,/ 
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Abstract .  A  concept  of  controllability  for  systems 

Ax  +  Bx  =  Cu(t)  in  which  A  may  be  singular  is  introduced. 

-  /A 

When  det(As+B)  0 ,  s  G  (E,  this  is  shown  to  be  equivalent 

to  the  condition  that  c  (As+B)  C  =  0  implies  c  =  0.  It 

is  also  shown  that  when  such  a  system  is  controllable  and  C 

is  a  column  vector,  then  there  exists  a  feedback  u  =  g^x  +  g^x 

T 

such  that  A  -  Cgj^  is  non-singular  and  all  solutions  of 

•  T  T  * 

Ax  +  Bx  =  C(gQX+g|x)  decay  exponentially. 


where  x  and  u  are  IR*'-  and  ■lR”’-valued  functions,  respec¬ 
tively,  of  t  e  IR  and  where  the  constant  matrices  A  and 


B  are  n  x  n  and  C  is  n  x  ,j  and  all  have  real  elements. 

We  shall  say  the  system  is  s ingular  when  A  is  singular  and 
regular  when 

(1.2)  A(s)  =  det(As+B)  0,  s  e  IR  (or  C) . 

Condition  (1.2)  is  the  criterion  for  the  pencil  of  matrices 
(As+B)  to  be  regular  in  the  terminology  of  [2],  p.25. 

Our  results  for  regular  singular  systems  (1.1)  reduce  to  well- 
known  facts  for  the  case  when  A  is  non-singular.  Regular  singular 
systems  of  the  form 

(1.3)  Ax  +  Bx  =  f(t) 

have  been  treated  elsewhere  [1]  where  a  formula  for  the  solutions 
of  (1.3)  is  given.  We  give  an  alternative  development  and  an 
equivalent  formula  in  §2.  A  more  detailed  exposition  is  in  [3]. 

It  is  known  and  will  become  evident  in  §2  that  when  a  solution 
of  (1.1)  or  (1.3)  exists  for  given  initial  conditions,  then  the 


2. 


solution  is  unique.  When  £  in  (1.3)  or  u  in  (1.1)  is  con¬ 
tinuous  on  an  interval  y,  then  a  solution  x  ovi  y  will  be 
understood  to  be  a  function  which  is  differentiable  on  y  and 
which  satisfies  the  equation  everywhere  on  y.  Thus,  in  treating 
the  concept  of  controllability  in  regard  to  (1.1),  we  restrict 
the  controls  u  to  the  class  of  continuous  lR"’-valued  functions 
and  require  differentiability  of  the  corresponding  responses  x. 

We  shall  see  that  when  control  can  be  effected,  the  function  u 
can,  in  fact,  be  chosen  in  a  class 

(1.4)  =  {u:[tQ,»)  -»-1R"':u^^^  is  continuous,  k=0,l,...,y} 

for  some  u  >  0. 

Definition  1.1. 

System  (1.1)  is  controllable  (at  time  t^)  if  for  every 
€1R”,  there  is  a  u  €  “STq  and  a  x  >  0  such  that  there  is 
a  solution  x  of  (1.1)  satisfying  x(tp)  =  ^  and  x(tQ+T)  = 
it  is  controllable  from  zero  if  the  same  is  true  with  the  restric¬ 
tion  ?  =  0. 

Remark  1.2. 

It  should  be  clear  that  (1.1)  is  controllable  at  time  t^ 
if  and  only  if  it  is  controllable  at  time  0.  Since  A  and  B 
are  constant,  one  need  only  translate  the  control  which  effects 
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the  transfer  from  C  to  Accordingly,  unless  otherwise  noted,  we 

take  tQ  =  0  in  (1.4)  and  Definition  1.1. 

Our  main  result,  proved  in  §3,  is  the  following. 


Theorem  1.3. 

Let  (1.1)  be  regular.  Then  (1.1)  is  controllable  if  and  only 
if 

(1.5)  c  €]R",  c^(As+B)'^C  =  0,  s  €1R,  implies  c  =  0. 

Moreover,  when  (1.1)  is  controllable,  we  may  take  x  >  0 
arbitrarily,  regardless  of  C,^. 

Remark  1.4. 

If  the  coefficient  matrices  A,B,C  have  elements  in  C,  the  complex 
numbers,  then  our  results  are  valid  also  in  the  corresponding 
context;  that  is,  we  may  replace  IR*'  and  ]r"*  in  (1.4),  Definition 
1.1  and  (1.5)  by  and  c"*,  respectively,  with  x  being 

C^-valued.  It  will  be  clear  that  our  proofs  remain  valid  in  this 
context  when  appropriate  trivial  modifications  are  made.  It  should 
be  noted  that  (As+B)’^  is  not  generally  defined  for  all  s  e  IR 
(or  C) .  However,  when  (1.2)  holds,  the  elements  of  (As+B)’^  are 
rational  functions  of  s  and  (1.5)  is  to  be  interpreted  to  mean 
that  the  elements  of  c  (As+B)  C  are  each  the  zero  function  at 


points  s  where  they  are  defined.  Note  that  it  is  irrelevant 


whether  s  is  interpreted  as  a  variable  over  IR  or  over  C. 

Remark  1,5. 

Singularity  of  A  is  not  an  hypothesis  in  Theorem  1,3. 

When  A  is  non-singular  it  is  easy  to  show  by  expanding  (As+B) 
in  powers  of  s  ^  for  large  s  that  (1.5)  is  equivalent  to  the 
condition 

rank[Bj^,A^B^,.  .  .  ,aJ'^B^]  =  n,  B^  =  A'^C,  A^  =  -A'^B. 

This  is  the  well-known  condition  (Theorem  5,  p.81,  [5])  for 
controllability  of  the  system  x  =  A^x  +  B^u  equivalent  to  (1.1) 
when  A'^  exists. 

In  §4  we  prove  a  stabilizability  result  for  (1.1).  It  is 
convenient  to  use  the  following. 

Definition  1.6. 

System  (1,1)  is  stabil izable  if  there  exist  real  m  x  n 
matrices  Gq  and  such  that  with 

(1.6)  u  =  Ggx  +  G^i 


in  (1.1)  the  resulting  system 
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(1.7)  Ax  +  Bx  =  0,  A  =  A  -  CG^,  B  =  B  -  CG^ , 

has  A  non-singular  and  all  solutions  tend  to  zero  exponentially 
as  t  ->■  +“,  or,  equivalently,  det(As+B)  has  degree  n  and 
all  its  zeros  have  negative  real. parts. 

The  main  result  in  §4  is  Theorem  4.1  which  states  that  if 

(1.1)  is  regular  and  controllable  and  m  =  1,  then  (1.1)  is 
stabilizable.  Corollary  4.2  affirms  the  analogous  result  for 
certain  cases  with  m  >  1.  We  intend  to  treat  the  general  case 
m  >  1  at  a  later  time. 

2 .  Decomposition  of  the  system. 

Throughout  this  section  we  assume  A  is  singular  and  we 
sketch  briefly  an  analysis  leading  to  an  explicit  formula  for 
the  solutions  of  (1.1)  in  this  case.  Additional  details  can  be 
found  in  [3].  Condition  (1.2)  implies  the  existence  of  a  Laurent 
expansion  for  (A+zB)’^  in  a  deleted  neighborhood  of  zero;  that 
is, 

CO 

(2.1)  (A+zB)'^  =  I  z'^Q.  ,  zee,  0  <  |z|  <  6 

k=-p  ^ 

for  some  6  >  0.  Here  1  since  A  is  singular  and 

the  n  X  n  matrices  Qj^  have  real  elements  when  A  and  B  do. 
(Many  of  the  relations  which  follow  appear  in  [4]  but  derived 
from  a  different  point  of  view.  It  should  be  noted  that  we  used 
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y  +  1  there  in  place  of  y  here.) 

From  (2.1)  it  readily  follows  that 

V  =  -Qk-l®’  =  -^Qk-l’  ^  ^  ° 

(2.2) 

QflA  *Q.iB  -  I„  .  AQ„  *  BQ_J  = 

where  1^^  denotes  the  n  ^  n  identity  matrix.  One  may  now 
show  that 

(2.3)  AQj^B  =  BQj^A,  k  >  -y, 

and  that 

k  <  -1 ,  j  ^  0 
k  >  0 ,  j  ^  0 
k£-l,  j£-l. 

From  these,  (2.2)  and  (2.3)  it  follows  that 
f2-4)  QgAQo  =  Qp  ,  Q.jBQ.j  =  Q_j 

and  since  Q_j^  f  0,  that 

(Q.jA)^  =  0  ,  (Q.^A)*"’^  i  0. 


(2.5) 
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If  we  define 


(2.6) 


then  from  (2.2)  and  (2.5)  we  see  that  Pq  and  are  com¬ 

plementary  projections; 


(2.7) 


0,1; 


Pq  *  ^1 


If  we  let  r  =  rank  Pg ,  then  r  <  n  since  A  is  assumed  to  be 
singular.  For  convenience  we  assume  r  >  0;  the  case  r  =  0  is 
included  in  what  follows  if  one  omits  various  terms  which  are 
vacuous,  in  effect,  in  that  case.  In  a  similar  way  the  case 
r  =  n  (A  non- singular)  is  included  in  what  follows.  Accordingly, 
we  define  p  =  n  -  r  and  consider  that 

(2.8)  r>0,  p=n-r>0. 


Now  let  X  be  n  x  r  and  Y  be  n  x  p  with  the  columns 
of  X  and  Y  forming  bases  for  the  ranges  (column  spaces)  of 
Pg  and  Pji  ,  respectively.  By  (2.7)  the  n  x  n  matrix 


(2.9) 


T  =  [X,Y1 


is  non-singular  and  we  define  r  x  n  and  p  x  n  matrices 
and  V,  respectively,  by 

(2.10)  T'^  = 


The  following  relations  then  hold: 

(2.11)  UX  =  UY  =  0,  VX  =  0,  VY  =  ,  XU  +  YV  =  1^^. 

Since  PqX  =  X  and  P^^Y  =  Y  relations  (2.11)  imply 

(2.12)  Pq  =  XU,  Pj  =  YV 

(2.13)  UPq  =  U,  UP^  =  0,  VPp  =  0,  VP^  =  V. 

Below  we  shall  need  the  r  x  r  matrix  b  ^nd  p  x  p 
a  defined  by 

(2.14)  B  =  -UQgBX  ,  a  =  -VQ_^AY. 


fi 


Using  (2.3)  and  (2.4),  one  finds  that  YaV  = 
then  implies 


a 


P  = 


.P-1 


0. 


matrix 

.5) 

. J 


(2.15) 


0 


f 
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We  shall  also  need  the  n  x  n  matrix 


(2.16) 


S  = 


UQr 


VQ 


-1 


Observe  that  [AX,BY]S  =  APqQq  +  ^^i^-l  =  (2-12),  (2.4) 

and  (2.2).  Hence 


(2.17) 


S'^  =  [AX, BY] . 


Lemma  2.1. 

Let  (1.1)  be  regular  and  suppose  (2.8)  holds.  Then  (1.1)  is 
equivalent  to 


(2.18p) 

(2. ISj) 


where 


Xq  -  6Xp  -  TpU 


^1  '  ‘^^l  ~  ^l'^ 


(2.19) 


Xxq  +  Yxj,  Xq  e  ir’^,  Xj^  e  ]rp 


and 


♦ 


(2.20) 


Fq  =  UQpC 


Fi  =  VQ.^C. 
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Proof : 

Let 


(2.21) 


T‘^x 


so  that 


(2.22) 


Xp  =  Ux,  X,  =  Vx. 


Using  the  substitution  (2.19),  we  see  that  (1.1)  is  equivalent  to 

(2.23)  (SAT)x  +  (SBT)x  =  SCu. 

Computing  the  submatrices  comprising  SAT  we  find 

UQqAX  =  I^,  UQqAY  =  0,  VQ_j^AX  =  0,  VQ_^AY  =  -a 

by  virtue  of  (2.3),  (2.11),  (2.12),  (2.13)  and  (2.14).  Similar 
computations  give  the  submatrices  comprising  SBT  and  we  have 


0 

"-B  0 

r 

,  SBT  = 

0 

-a 

(2.24) 


SAT 


1 1 . 

It  is  now  oviilont  that  (1.1)  is  etjiiivalont  to  (2. IS.),  i  =  0,1, 
l>y  virtue  of  (J.21)  aiul  [2.20).  (#J 

Lemma  2.2. 

I.et  Xj  lie  a  solution  of  {2.IS|)  on  an  interv.il  ^  with 
u  eontinuous  on  f.  If  we  ilel'ine 

(2.  2.S)  W|^  =  a''  ’'xj  ,  k  =  0,1 . II. 

then  W|^  exists  on  ^  aiu! 

1 1 "  k. '  I  * 

(2.2(1)  Wj,^j  =  a'  I'lU  +  W|^,  K  -  1>,...,))-1, 
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which,  in  effect,  is  (2.26).  Since  Wp  =  =  0  by  (2.15),  we 

see  that  Wp  is  uniquely  determined  by  u.  The  same  follows  in¬ 
ductively  for  Wj^,...,w^  by  (2.26).  With  k  =  0  in  (2.26)  we 
have  Wj  =  ^^1^  which  is  (2.27)  for  j  =  1.  An  inductive 

proof  using  (2.26)  then  establishes  (2.27)  for  any  j,  1  £  j  £  P , 
when  u  is  j  -  1  times  differentiable.  ^ 

Remark  2.3. 

Taking  k  =  u  in  (2.25),  we  see  that  the  previous  lemma 
implies  that  is  a  formula  for  the  solution  of  (2.18j^) 

when  u  is  continuous.  For  example,  if  p  =  3  one  gets  from 
(2.26)  and  the  fact  that  Wp  =  0,  that 

Xi  =  W3  =  Fju  +  ^(ar^u+^^r^u) 


with  no  implication  that  u  and  u  exist.  When  u  is  y  -  1 
times  differentiable  on  then  (2.27)  with  j  =  y  gives  the 

easily  written  formula 

(2.28)  X,  (t)  =  a^r,u^^^(t),  t  e  I. 

^  k=0  ^ 

In  any  case,  with  u  continuous  on  _f  the  solution  of  (2.18p) 
with  Xp(0)  given  is 
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(2.29) 


XgCt) 


e(t-s) 


rQU(s)ds . 


Observe  that  whereas  Xq(0)  can  be  chosen  arbitrarily  and  in¬ 
dependently  of  u,  the  value  of  Xj(0)  is  determined  by  u. 
Formulas  (2.28)  and  (2.29)  are  equivalent  to  equation  (19),  p.419, 
in  [1]. 


3 .  Controllability  of  the  system. 

Controllability  for  (1.1)  clearly  implies  that  (2.18^)  and 
(2.18j^)  are  each  controllable.  Whereas  criteria  for  controllability 
of  (2.18p)  are  well-known,  this  seems  not  to  be  the  case  for  (2.18^^). 
In  connection  with  the  latter  the  following  is  important. 

Lemma  5.1.  Let  be  a  solution  of  (2.18^)  on  an  interval  y 

where  u  is  continuous.  If  Cj  6  IR*^  satisfies 

(3.1)  ~  ^  “  0,...,y-l, 

then 


(3.2) 


CjXj(t)  =  0 


on  y. 
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Proof; 

When  u  is  y  -  1  times  differentiable,  this  follows 
immediately  from  (2.28).  If  u  is  merely  continuous,  then  (3.2) 
follows  from  (3.1)  by  an  induction  argument  employing  (2.26)  in- 
as  much  as  =  0  and  w^  = 


Theorem  3.2. 

Let  (1.1)  be  regular.  In  order  that  it  be  controllable  it 
is  necessary  that 

(3.3)  rank  Hq  =  r,  rank  =  p 

where 

(3.4)  Hq  =  [ro,erg,...,B’"irQ],  ii^  =  [r^,ar^ . 


Proof : 

It  was  pointed  out  above  that  if  (1.1)  is  controllable,  then 
both  (2.18q)  and  (2.18j^)  are  controllable.  We  must  then  have 
rank  Hq  =  r.  The  condition  rank  llj  =  p  follows  from  Lemma  3.1.  Indeed  if 
rank  <  p  then  there  is  a  C  IR^,  Cj^  0  such  that  =  0. 

But  then  by  (3.2)  we  have  Cj^Xj(t)  =  0  so  there  can  be  no  control 
u  €*2^0  transferring  x  from  x(0)  =0  to  x(t)  =  Yc^^  for  t  >  0 
inasmuch  as  then  Xj(t)  =  Cj  (cf.  (2.19)  and  (2.21))  for  which 
c{xi(T)  ^0.  g 
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Theorem  3.3. 

Let  (1.1)  be  regular.  Then  it  is  controllable  from  zero  if 
and  only  if  (3.3)  holds.  Moreover,  when  (3.3)  holds,  then 
T  >  0  in  Definition  1.1  can  be  chosen  arbitrarily  and  one  can 
choose  u  €  (tQ=0)  such  that  . 

(3.5)  u(0)  =  u'(0)  =  =  u^‘^^(0)  =  0 

so  that 

(3.6)  x(0)  =  0  when  x(0)  =  0. 

Proof : 

The  proof  given  for  Theorem  3.2  suffices  to  show  that  (3.3) 
is  necessary  for  (1.1)  to  be  controllable  from  zero.  To  prove 
that  (3.3)  is  sufficient  we  restrict  ourselves  to  controls 
u  €  defined  by 

ft 

(3.7)  u(t)  =  (t-s)^v(s)ds ,  t  >  0 

Jo 

where  v  €  ‘2rg(tQ=0).  For  such  u  we  have 

(3.8)  u^^^(t)  =  (p)j^  (t-s)^'*^v(s)ds ,  k  =  0,1,. ..,n 

where  (p)j^  =  y:/(M-k):  so  that  (3.5)  holds. 
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For  a  solution  of  the  equivalent  system  (2.18^),  i  =  0,1, 
we  then  necessarily  have  Xj^(O)  =  Xj(0)  =  0  by  (2.28).  When 
x(0)  =  0,  then  Xg(0)  =  Ux(0)  =0  so  by  (2.29)  and  (3.7)  we 
have 


(3.9) 


XgCt) 


•t 

^  Kg(t-s)rgV(s)ds,  t  >  0 


where 


(3.10) 


Kg(t)  =  e^*^^'®^a^da 


t  >  0. 


Observe  that  then  Xg(0)  =  0.  Hence  for  controls  as  in  (3.7)  we 
have  (3.6). 

Substituting  (3.8)  into  (2.28),  we  may  write 

(3.11)  Xj^(t)  =  I  Kj  (t-s)  rj^v(s)ds ,  t  >  0 
where 

(3.12)  K,(t)  =  (u),t'^’V. 

^  k=0 


Combining  (3.9)  and  (3.11)  and  using  (2.19),  we  have 
(3.13)  x(t)  =  T  f  K(t-s) rv(s)ds 

Jo 
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with 


o 

o 

1 _ 

(3.14) 

K(t)  = 

> 

r  = 

r 

0  K^(t) 

L  ij 

We  now  define  a  n  x  n  matrix  function  W  by 


(3.15) 


W(T) 


K(T-s)rr^K^(T-s)ds. 

■’0 


It  is  evident  that  W(t)  is  real  symmetric  ajid  positive  definite  or  semi- 
definite.  Moreover,  W(t)  is  singular  for  x  >  0  if  and  only 
if  there  is  some  c  6  ]R^,  c  0 ,  such  that 

(3.16)  c^KCt)r  =0,  0  £  t  <  T. 


If  W(t)  is  non-singular  for  some  x  >  0,  then  for  any  ?  e  ]R” 
we  can  define  v  by 


(3.17) 


v(t)  =  r'^K^(x-t)W'^(T)T'^C. 


With  this  V  in  (3.7)  the  resulting  control  u  G  »  satisfies 
(3.5)  and,  by  virtue  of  (3.13),  gives  x(x)  =  c  when  x(0)  =  0 
as  well  as  x(0)  =  0.  The  theorem  will  be  proved  then  when  we 
show  that  (3.3)  implies  that  W(x)  is  non-singular  for  all  x  >  0. 
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This  is  included  in  the  following  lemma. 

Lemma  3.4. 

For  any  x  >  0,  W(t)  defined  in  (3.15)  is  non-singular  if 
and  only  if  (3.3)  holds. 

Proof : 

Let  T  >  0  be  given  and  let  c  £  satisfy  c^W(t)  =  0. 
This  is  equivalent  to  (3.16)  which,  in  turn,  is  equivalent  to 

(3.18)  c^K(*^^(0)r  =  0  ,  k  =  0,1,2,... 


inasmuch  as  K(t)  is  an  entire  function.  From  (3.10)  one  can 
compute  that 


K 


(k) 

0 


(0) 


fo. 


k  =  0,1, ...  ,y 
k  >  y  +  1 . 


From  (3.12)  we  get 


K 


(k) 

1 


(0) 


y !  a 


y-k 

> 


0, 


k  -  l,...,y 
L  =  0  or  k>y-*-l. 


T  T  T 

Accordingly,  if  we  take  c  =  [c^,  Cj] 
then  (3.18)  is  equivalent  to 


where  Cq  €  ,  Cj^  £  IR^, 
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(3.19) 


cJe^Fo  =  0 

T  k 

cja  Ti  =  0 


k  =  0,1,2,... 
k  =  0,l,...,y-l 


These  in  turn  are  equivalent  to 


(3.20) 


T 

c{Hi 


0 


T 

where  Hq  and  are  defined  in  (3.4).  If  c  W(t)  =  0  and 

c  0 ,  then  (3.20)  holds  with  Cq  0  or  Cj  ^  0  so  (3.3)  fails. 
Conversely,  if  (3.3)  fails,  then  (3.20)  holds  for  some  Cq,Cj  not 
both  zero.  The  corresponding  c  0  and  satisfies  c  W(t)  =  0. 
This  proves  the  lemma  and  completes  the  proof  of  Theorem  3.3.  ^ 

Lemma  3.5. 

Conditions  (1,5),  (3.3)  and 

(3.21)  c  €  IR*^,  c^(-As+B)'^C  =  0,  s  e  IR,  implies  c  =  0, 
are  all  equivalent. 

Proof : 

The  equivalence  of  (1.5)  and  (3.21)  is  obvious  since  -s 
may  be  substituted  for  s  in  either.  To  show  that  (1.5)  and  (3.3) 
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are  equivalent  we  use  (2.24)  to  note  that 


(3.22)  S(As+B)T  =  diag (sl^- B  ,  I p- sa) . 


It  follows  that 


(3.23)  J(As+B)'^C  =  cJ(sI^-B)'^rQ  +  c]^(Ip-sa) 
where 

C3.24)  cj  =  Jx.  c}  =  Jy,  J  =  [cJ.cJlI-1. 


But  for  large  |sl,  s  €  IR,  we  have 


(si  -6)'^  =  I  b’ 

^  k=0 


while  for  any  s  €  IR  we  have 


(I  -sa)'^  =  ^'l 
P  k=0 


by  virtue  of  (2.15).  From  (3.23)  we  see  then  that  c^(As+B)'^ 
s  €  IR,  is  equivalent  to  (3.19)  and  hence  to  (3.20).  By  virtu 

(3.24)  it  is  evident  then  that  (1.5)  holds  if  and  only  if  (3.3 
holds.  This  completes  the  proof  of  the  lemma.  @ 


C  =  0, 
e  of 


M 

J 
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We  are  now  in  a  position  to  prove  our  main  result: 

Proof  of  Theorem  1.5: 

As  was  pointed  out  earlier  we  may  assume  A  is  singular 
and  (2.8)  holds.  Hence  if  (1.1)  is  regular  and  controllable, 
(3.3)  must  hold  by  Theorem  3.2  and  this  implies  (1.5)  by  Lemma 
3.5.  Now  suppose  (1.5)  holds.  Then  (3.3)  holds  so  by  Theorem 
3.3  system  (1.1)  is  controllable  from  zero  with  any  x  >  0  and 
with  u  £  satisfying  (3.5)  and  the  resulting  solution  x 
satisfying  (3.6).  But  (1.1)  also  implies  (3.21)  by  Lemma  3.5 
so,  by  the  same  argument  just  given,  the  system 

(3.25)  -Ax  +  Bx  =  Cu 

is  controllable  from  zero  with  any  x  >  0  and  analogous  u  € 
and  solution  x.  Hence,  if  (1.1)  holds  and  GIR’^,  x  >  0  are 

given,  there  is  a  control  u^  €  steering  the  corresponding 

solution  x^  of  (3.25)  from  x^(0)  =0  to  x^(x/2)  =  C  and 
there  is  a  control  u  6^^  steering  the  corresponding  solution 
x^  of  (1.1)  from  x^(0)  =0  to  x^(t/2)  =  C.  The  controls  u^ 
and  u^  satisfy  (3.5)  and  by  (3.6)  we  have  x^(0)  =  0,  x^(0)  =  0. 
If  we  define 

f  u^(T/2-t)  ,  0  <  t  <  T/2 

u(t)  =  <;  2 

[  u^(t-T/2) ,  t  1  f/2. 


(3.26) 
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then  u  £‘2;'^.  The  function  x  defined  by 


C3.27) 


x(t) 


X^(T/2-t) 
x^  (t-T/2) 


0  <  t  <  t/2 
t  >  t/2 


is  continuous  and  differentiable  for  t  >  0  (even  at  t  =  t/2) 
and  is  a  solution  of  (1.1)  with  u  given  by  (3.26)  which  satisfies 
x(0)  =  C  and  x(t)  =  z,.  This  completes  the  proof.  ^ 

4 .  Stabilizability  of  the  system. 

When  (1.1)  is  regular,  the  solutions  of  the  homogeneous 
system 


(4.1) 


Ax  +  Bx  =  0 


can  be  found  by  setting  u(t)  so  in  (2.28)  and  (2.29).  One  gets 
XQ(t)  =  e^^Xg(O),  Xj(t)  s  0,  and  by  (2.19),  we  find  the  solutions 
of  (4.1)  in  the  form 


(4.2) 


x(t)  =  Xe^^Cg  ,  Cq  ^ 


Thus  the  initial  values  are  restricted  to  the  range  of  P,,  by 


x(0)  = 


0 


(4.3) 
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inasmuch  as  the  columns  of  X  are  a  basis  for  the  range  of  Pg. 

By  (2.15)  we  have  det(Ip-sa)  =1  so  from  (3.22)  we  see 
that  the  eigenvalues  of  B  are  precisely  those  X  for  which 
A(X)  =  0  (cf.  (1.2)).  The  polynomial  A  is  of  degree  r  and 
the  solutions  (4.2)  form  an  r-dimensional  space  over  IR.  If  all 
eigenvalues  X  have  negative  real  parts,  then  the  zero  solution 
of  (4.1)  is  asymptotically  stable  relative  to  the  allowable 
initial  values  x(0)  given  in  (4.3).  However,  a  bounded  input 
u  in  the  system  (1.1)  need  not  result  in  a  bounded  response  x. 
Indeed,  if  r  <  n,  u  ^  2  and  u  €  ,  we  see  from  (2.28)  that 

then  Xj^  and  hence  x  may  be  unbounded  when  is  un¬ 

bounded  . 

In  regard  to  (1.1)  we  ask  then;  what  conditions  assure  that 
the  system  can  be  stabilized  by  a  linear  feedback  control  so  that 
the  zero  solution  of  the  combined  system  plus  feedback  is 
asymptotically  stable  relative  to  arbitrary  initial  conditions 
in  the  state  space  IR*^?  It  should  be  clear  that  when  r  <  n,  a 
feedback  u  =  Gx  with  G  being  n  x  m  will  not  suffice.  The 
combined  system  is  then  Ax  +  (B-CG)x  =  0  and  the  degree  of 
det(As+B-CG)  will  still  be  less  than  n  since  A  is  singular. 
The  corresponding  solutions  and  allowable  initial  conditions  will 
be  constrained  to  lie  in  a  proper  subspace  of  IR*^.  To  achieve 
stabilization  when  A  is  singular  the  feedback  must  contain  a 
term  involving  x.  Thus  we  ask  for  conditions  which  assure 
stabilizability  in  the  sense  of  Definition  1.6.  Our  principal 
result  in  this  direction  is  the  following  theorem. 
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Theorem  4.1. 

Suppose  the  system  (1.1)  is  regular  and  that  C  is  n  x  1 
(m  =  1).  If  (1.1)  is  controllable,  then  it  is  stabilizable. 

Proof : 

As  before,  we  treat  (1.1)  in  the  decomposed  form  (2.18^), 
i  =  0,1,  under  the  assumption  (2.8);  the  extreme  cases  r  =  0 
and  r  =  n  (A  non-singular)  are  included  thereby  when  vacuous 
terms  are  omitted.  By  hypothesis  (1.1)  is  controllable,  so  (1.5) 
holds  by  Theorem  1.3  and  hence  (3.3)  holds  by  Lemma  3.5.  Since 
Tq  is  r  X  1  (m  =  1) ,  then  Hg  is  r  x  r  and  non-singular. 

It  follows  (Theorem  7,  p.90,  [5])  that  there  exists  a  non-singular 
r  X  r  matrix  Fg  such  that  if 

(4.4)  yg  =  FgXg, 

then  (2.18g)  is  equivalent  to 

(4.5)  Xg  -  M(3)yg  =  e^u 

where  is  the  r^  column  of  and  M(B)  is  a  companion 

matrix  for  the  characteristic  polynomial  of  8;  specifically. 


I  1 . 0 1 


M{|0 


{)  1  0 

I 

(1  0  0 


0 

0 

1 


who  VO  vU't(.  \lj. 

Tho  ni.itrix  II  j 
siiioo  r;ink  II ^  p, 
wo  IIUISI  li.ivo  |l  •  (> 

sino.iilar.  Thoro  is 
t  ha  t  i  r 


.\  ,  A  . 

V  I 

i  II  (  .A ,  I )  i  s  p  ''  |i 
Siiuo  It  is  p  ^  p 
lloiuo  II  ‘  p  aiul 
I  ho  11  a  non  '  s  i  nt'.n  I  a  V 


so  wo  mipit  havo  p  ■  p 
aiul  sal  i  s  r  ios  l .  I  S  ) 

II I  is  sinKiio  aiul  noil 
p  s  ,>  matrix  1^  siioli 


I  •» .  7  ) 


I|X| 


thou  (J.ISj)  is  otm  i  va I  on  I  to 

(  1  .  S )  ♦  Vj  -  0^11 


wlioro  o  is  tlio  ('**’  ooliniin  ol  I  aiul  M|it)  a  oonipanion 
P  p 

matrix  lor  (ho  olia  rao  t  or  i  s  (  i  o  polynomial  of  u  ;  spoo  i  T  i  i  a  I  I  y  , 


Milt)  is  oT  tho  samo  rorm  as  M(lO 
row  is  7.0 ro.  I.ottiny. 


in  (l.li)  oxoopt  that  tho  last 
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(4.9) 


and  combining  (4.5)  and  (4.8),  we  have  the  system 


Pr  ® 

-M(0)  0 

e 

(4.10) 

y  + 

y  = 

T 

0  -  M(a) 

0  I 

e 

L  p_ 

which  is  equivalent  to  (1.1)  under  the  hypotheses  of  the  theorem. 

It  is  convenient  now  to  replace  (4.10)  by  the  system  equiva¬ 
lent  to  it  obtained  by  subtracting  the  last  equation  of  the  system 
from  the  r^*^  equation.  Since  the  last  row  of  M(ct)  is  zero,  this 
is 


(4.11) 


0 

M(a) 


'M(B)  -  Dp 

0  ' 

y  + 

y  = 

0  I 

e 

L  pj 

L  pJ 

where  is  r  p  and  all  its  elements  are  zero  except  for  a 

one  in  the  r^^  row  and  column. 

T 

Next  let  ej^  =  [1,0,. ..,0]  be  I  x  p  and  replace  u  in 
(4.11)  by 


(4.12) 


-  v. 


The  result  may  be  written,  after  changing  signs  in  the  last  p 


I 
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equations  of  the  system,  in  the  form 


We  now  define  two  permutation  matrices  and  :  tTq  is 

p  p  and  reverses  the  order  of  the  columns  of  a  matrix  (of  p 
columns)  when  used  as  a  postfactor;  tTj  is  p  P  and  reverses 
the  order  of  the  first  p  -  1  rows  and  leaves  the  last  row  un¬ 
affected  when  used  as  a  prefactor  on  a  matrix  with  p  rows.  Note 
2  2 

that  11^  =  17^^  =  Ip.  In  (4.13)  we  make  the  substitution 

(4.14)  y  =  [diag(I^,Trg)  ]  z 


and  multiply  the  result  by  the  matrix  diag(I^,TTj^)  .  Since 

T 

TTj^ep  =  ep  and  (M(a)+epe2^)  =  Ip,  the  resulting  equivalent 
system  takes  the  form 

(4.15)  z  =ja^z  +  bv 


where 

(4.16) 


28. 


Note  that  all  elements  of  DqTTq  are  zero  except  for  a  one 

in  the  r^^  row  and  first  column..  Also,  the  element  in  the  i^^ 
t  h 

row  and  i  +  1  column  of  i  =  l,...,p-l,  is  a  one  and  all  others 

are  zeros  except  for  a  one  in  the  p  row  and  first  column. 

Hence,  from  the  form  of  M(8)  in  (4.6),  we  see  that  has  ones 
just  above  the  main  diagonal  and  all  zeroes  everywhere  above  these. 
It  is  then  clear  from  the  form  of  b  in  (4.16)  that 

(4.17)  H  =  [b,j^,...  ^'h] 

is  non-singular.  Hence  there  exists  a  feedback  (Theorem  9,  p.97, 

[5]) 

(4.18)  V  =  g'^z  ,  g  €  IR^ 

which  when  substituted  into  (4.15)  gives  the  system 

(4.19)  z  =  Cc/+bg’^)z 

with  Si/  +  bg*  stable;  that  is,  g  can  be  chosen  so  that  all 
eigenvalues  of  -Si/  +  bg^  have  negative  real  parts  and  all  solutions 
of  (4.19)  tend  to  zero  exponentially  as  t  +“. 

The  substitutions  (4.4)  and  (4.7)  can  be  written  as 

y  =  [diag(FQ,Fj)]T'^x 


(4.20) 
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and  (4.14)  can  be  written 


(4.21) 


z  =  [diag(I^,TTQ)  ]y 


since  =  Ip.  The  combined  feedback  resulting  from  (4.12)  and 
(4.18)  is  thus  of  the  form 


(4.22) 


T-  T  T  T* 

u  =  -  g  z  =  gpx  +  g^x 


for  some  gg*®!  specifically. 


(4.23) 


J  -  T 
§0  “  "8 


.^FiV 


T  T 


Using  (4.20)  and  (4.21),  we  can  express  (4.19)  in  terms  of  x. 

The  resulting  equation,  when  multiplied  by  the  inverses  (in  reverse 
order)  of  the  several  matrices  used,  in  effect,  as  prefactors  in 
going  from  (1.1)  to  (4.19),  will  produce  the  form  (1.7)  with 
Gq  =  gj,  Gi  *  g]^;  that  is,  (1.7)  with 

(4.24)  A  =  A  -  Cg]^  ,  B  =  B  -  CgJ. 

The  term  Ax  arises  from  the  term  z  =  I^z  in  (4.19)  so  A  must 

n 

be  non-singular  and  all  solutions  x(t)  of  (1.7)  tend  to  zero 
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exponentially  as  t  ->  »  since  all  solutions  z(t)  of  (4.19) 
do  so.  This  completes  the  proof,  of  the  theorem.  0 

In  the  above,  C  was  assumed  to  be  n  x  1.  We  can  use  this 
to  prove  the  stabilizability  of  (1.1)  for  some  cases  when  C  is 
n  X  m  with  m  >  1. 

Corollary  4.2. 

Suppose  the  system  (1.1)  is  regular,  that  rank  A  =  n  -  1 
and  that  the  zeros  of  A(s)  =  det(As+B)  are  distinct  and  none  is 
zero.  If  (1.1)  is  controllable,  then  it  is  stabilizable. 

Proof ; 

Since  rank  A  =  n  -  1,  A  is  singular  so  the  polynomial 
A(s)  has  degree  r  <  n.  (A(s)  ^  0  since  (1.1)  is  regular.)  We 
may  write 

(4.25)  A(s)  =  dgS^+.-.+d^,  dp  0,  d^  0, 

the  condition  d^  0  being  a  result  of  the  hypotehsis  that 
A(0)  ^  0.  From  (4.25)  we  find  (z  £  C) 

det(A+zB)  =  z"a(1/z)  =  z"‘^ (d^z’^+ .  .  .  +  dp) . 

The  cofactor  of  at  least  one  element  in  A  +  zB  is  non-zero 
at  z  =  0  since  rank  A  =  n  -  1.  Hence  since  dp  0 


31. 


lim  z’^‘’’(A+zB)‘^ 
z->0 


exists  and  is  non-zero  so  (A+zB)'^  has  a  pole  of  order  n  -  r 
at  z  =  0.  Thus  y  in  (2.1)  satisfies  y  =  n  -  r  =  p.  Since 
r  <  n,  we  have  p  >  0  and,  again  for  notational  convenience 
only,  we  assume  r  >  0  so  that  (2.8)  holds. 

Now  let 


(4.26) 


0 

a 


where  3  and  a  are  defined  in  (2.14)  as  in  the  development  in 
§2.  Consider 


(4.27) 


H  =  [r,Er, . . . ,E"'ir] 


where  r  is  as  in  (3.14).  Since 


y  =  p  and  =  0 ,  we  find 


H  = 


»o  eX 
0 


where  =  [Tp  jBTp , .  .  .  ,B^  ^^0^  ^*1 

(3.4).  But  the  eigenvalues  of  g  are  the  zeros  of  A(s)  so  our 

hypotehses  imply  that  g  is  non-singular.  The  controllability 
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of  (1.1)  implies  (3.3)  by  virtue  of  Theorem  1.3  and  Lemma  3.5. 

It  follows  then  that  rank  H  =  n. 

Since  y  =  p  and  a  is  p  p  and  nilpotent  of  order 
y ,  the  Jordan  canonical  form  for  a  is  M(a)  as  described 
just  after  (4.8).  Any  Jordan  form  for  8  is  matrix  diagonal 
with  distinct  diagonal  elements  and  none  zero  by  virtue  of  our 
hypothesis  regarding  the  zeros  of  A(s).  Hence  in  any  Jordan 
form  for  E  no  two  Jordan  blocks  have  the  same  eigenvalue.  It 
follows  then  from  Theorem  6,  p.  86  of  [5]  that  there  exists  a 
c  €  IR™  such  that 

H  =  [r,Er,. . .  ,  r  =  rc 

has  rank  n.  This  in  turn  implies 

(4.28)  rank  Hq  =  r,  rank  =  p 

where  and  are  as  in  (3.4)  with  Tq  and  replaced 

by  TgC  and  T^c,  respectively. 

Again  using  Lemma  3.5  and  Theorem  1.3,  we  may  now  conclude 

that 

(4.29)  Ax  +  Bx  =  Cu  ,  C  =  Cc ,  u  €  ]R^ 

is  controllable.  Since  now  C  is  n  x  i ,  Theorem  4.1  implies 
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that  (4.29)  is  stabilizable .  But  (4.29)  is  (1.1)  with  u  =  cu 
and  the  stabilizing  feedback  u  -  g^x  +  gjx  for  (4.29)  determines 
a  stabilizing  feedback  (1.6)  for  (1.1)  with  =  cgj,  Gj^  =  cg^.  ^ 

It  may  be  noted  that  the  hypotheses  in  Corollary  4.2  do  not 
imply  r  =  n  -  1.  The  following  is  a  case  in  which  r  <  n  -  1. 


1 

1 

0 

0 

1 

5" 

A  = 

0 

1 

0 

> 

B  = 

0 

0 

1 

_0 

0 

0_ 

_1 

0 

0_ 

Here  rank  A  =  2  =  n  -  1  but  A(s)  =  s  +  1  has  degree  1  <  n  -  1. 
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